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Abstract 

In e + e~ collisions using the CLEO detector we have studied the decay of 
the D° to the final state K^tt + tt~ with the initial flavor of the D° tagged 
in charged D* decay. We use the Dalitz technique to measure the resonant 
substructure in this final state and clearly observe ten different contributions 
by fitting for their amplitudes and relative phases. We observe a K* + n~ 
component which arises from doubly Cabibbo suppressed decays or D — D° 
mixing. 
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Weak hadronic decays of charmed mesons are expected to proceed dominantly by resonant 
two-body decays in several theoretical models [1-5]. A clearer understanding of final state 
interactions in exclusive weak decays is an important ingredient for our ability to model 
decay rates as well as for our understanding of interesting phenomena such as mixing [6]. In 
this context an interesting final state is D° — > fr°7T + 7r~ which can proceed through a number 
of two-body states. Previous investigations [7-11] of the substructure in this channel were 
limited by statistics to the Cabibbo favored decays. A key motivation is to observe one or 
more of the D° — > K 7i + 7T^ resonant submodes that proceed via mixing or double Cabibbo 
suppression, such as K* + n~ or _K"o(1430) + 7r _ , and to measure their phase relative to the 
corresponding unsuppressed K°7r + 7r~ submodes. 

This analysis uses an integrated luminosity of 9.0 ftr 1 of e+e~ collisions at y/s 10 GeV 
provided by the Cornell Electron Storage Ring (CESR). The data were taken with the 
CLEO II. V configuration of the CLEO II multipurpose detector [12]. A silicon vertex detector 
(SVX) was installed in the upgraded configuration [13]. 

The event selection is similar to that used in our search for D°-D° mixing via the process 
D° — > D° — > K + n~ [14]. We reconstruct candidates for the decay sequence D* + — > 7T+.D , 
D° — > Ksit + it~ . Consideration of charge conjugated modes is implied throughout this Letter. 
The charge of the slow pion (nf or n~) identifies the charm state at t = as either D° or D°. 
We require the D* + momentum p^* to exceed 2.0 GeV/c, and we require the D° to produce 
the final state K s n + 7i~ . We reconstruct K s ^7r + 7r~ with the requirement that the daughter 
pion tracks form a common vertex, in three dimensions, with a confidence level > 10~ 6 . 
Signal candidates pass the vertex requirement with 96% relative efficiency Throughout this 
Letter, relative efficiency is defined as the number of events passing all requirements relative 
to the number of events when only the requirement under study is relaxed. 

Our silicon vertex detector provides precise measurement of the charged tracks in three 
dimensions [15]. We exploit the precision tracking of the SVX by refitting the Kg and ^ 
tracks with a requirement that they form a common vertex in three dimensions. We use the 
trajectory of the K s 7r + 7i~ system and the position of the CESR luminous region to obtain 
the D° production point. We refit the 7r+ track with a requirement that the trajectory 
intersect the D° production point. We require that the confidence level of each refit exceed 
10~ 4 . The signal candidates pass the D° mass and decay vertex requirement with 85% and 
91% relative efficiency, respectively. 

We reconstruct the energy released in the D* + ^7rfD° decay as Q = M*—M—m n , where 
M* is the reconstructed mass of the 7r^K s 7[ + 'K~ system, M is the reconstructed mass of 
the K s 7t + -K~ system, and m n is the charged pion mass. The addition of the D° production 
point to the 7r+ trajectory yields the resolution oq = 220±4 keV, where Oq is the core 
value from a fit to a bifurcated student's t distribution. We obtain a resolution on M of 
a M = 4.8 ±0.1 MeV and a resolution on m K o of a m = 2.4 ± 0.1 MeV. 

S K g 

We apply a set of 'prophylactic' requirements to exclude candidates with a poorly de- 
termined Q, M, or D° flight time, and Kg candidates that are likely to be background. 
The typical computed a for Q, M and flight time error is 150 keV, 5 MeV and 0.5 t d o, 
respectively. These are computed from the reconstruction covariance matrix of the daugh- 
ters of the D° candidate. We reject candidates where <tq, <tm or flight time error exceeds 
300 keV, 10 MeV or 2.0tdo, respectively. The relative efficiencies for the signal candidates 
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FIG. 1. Distribution of a) Q and b) M for the process D° — > i^7r + 7r . The candidates pass 
all selection criteria discussed in the text. 

to pass these cuts is 98%, 98% and 89%, respectively. We exclude Kg candidates with a 
vertical flight distance less then 500 microns to remove combinatoric background with zero 
lifetime. The signal candidates survive this requirement with 95% relative efficiency. The 
distributions of Q and M for our data are shown in Figure 1. 

We select 5299 candidates within three standard deviations of the expected Q, M, and 
m K o. The efficiency for the selection described above is nearly uniform across the Dalitz dis- 
tribution. In our simulation we generate the D° — > K s tc + tt~ uniformly populating the 
allowed phase space. We study our efficiency with a GEANT [16] based simulation of 
e + e~ — > cc events in our detector with a luminosity corresponding to more than three 
times our data sample. We observe that our selection introduces distortions due to inef- 
ficiencies near the edge of phase space, and fit the efficiency to a two dimensional cubic 
polynomial in [m 2 K0 _,rat + ) requiring that the efficiency does not change if the ir + and 7r~ 
are interchanged. Our standard result uses this efficiency parameterization to interpret the 
Dalitz distribution. To take into account a systematic uncertainty in our selection efficiency 
we compare the standard result with an efficiency that is uniform across the allowed Dalitz 
distribution. 

Figure 1 shows that the background is small, but non-negligible. Fitting the M dis- 
tribution to a signal shape as described above plus a quadratic background shape we find 
a background fraction of 2.1 ± 1.5%. We use this fraction as a constraint when fitting 
the Dalitz distribution. To model the background contribution in the Dalitz distribution 
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we consider those events in the data that are in sidebands five to ten standard deviations 
from the signal in Q and M and within three in m K o. There are 445 candidates in this 
selection, about four times the amount of background we estimate from the signal region. 
We compare these with background from our simulation also including e + e" annihilations 
producing the lighter quarks. We note that the background from the simulation is dom- 
inated by random combinations of unrelated tracks, and the shape of the background in 
the simulation and the sideband data sample agree well. The simulation predicts that the 
background uniformly populates the allowed phase space, and we model this contribution 
to the Dalitz distribution by fitting the data sideband sample to a two dimensional cubic 
polynomial in (m| r o 7r _,m^. +7r _). All parameters except the constant are consistent with zero. 
Other possible contributions to the background, resonances combined with random tracks 
and real D° decays combined with random soft pions of the wrong charge, are negligible in 
the simulation. The latter, called mistags, are especially dangerous to our search for "wrong 
sign" D° decays. Mistags populate the Dalitz distribution in a known way that depends 
on the shape of signal. When we analyze the Dalitz distribution we allow a mistag fraction 
with an unconstrained contribution. We have looked for the contribution of an anomalous 
resonance, such as p° or K*~ , plus random tracks to the background in the data by studying 
the sidebands in Q, M, and m K o, and conclude that any such contributions are negligible. 

Figure 2 shows the Dalitz distribution for the D° — > Kg-K + 7T^ candidates. A rich structure 
is evident. Contributions from K*~tt + and K° s p° are apparent. Depopulated regions exist 
suggesting destructive interference between some resonances and the dominant decay modes. 

We parameterize the Kgn + 7i^ Dalitz distribution following the methodology de- 
scribed in Ref. [17] using the same sign convention used in previous investigations 
of this decay channel [8,11]. We consider nineteen resonant subcomponents, K*~tt + , 
«;(800)-7r+, K*(UW)-ir + , K&1430)-tt + , ^(1430)-^+, K*(lQ80)-ir + , K 3 *(1780)-vr + , K° s p, 
K° s co, K° s p(U50), K°p(1700), A>(500), ATg/ (980), i^/ 2 (1270), i^/ (1370), i^/ (1500), 
1^/0(1710), and the "wrong sign" K* + n~ and i^Q(1430) + 7r~, as well as a non-resonant con- 
tribution. The parameters of the resonances are taken from Ref. [18] except for the /o's. We 
use Ref. [19] for the /o(980) and the coupled channel analysis of Ref. [20] for the /o(1370), 
/o(1500) and /o(1710). We consider that each of the resonances has its own amplitude and 
relative phase. The non-resonant contribution is modeled as a uniform distribution across 
the allowed phase space with a fixed relative phase. The phases and widths of the resonance 
contributions vary as given by the spin of the resonance as described in Ref. [17]. 

This study is sensitive only to relative phases and amplitudes. Thus we fix one phase 
and one amplitude. To minimize correlated errors on the phases and amplitudes we choose 
the largest color suppressed mode, K^p, which should be out of phase with the unsuppressed 
modes to have a fixed zero phase and an amplitude of one. Since the choice of normaliza- 
tion, phase convention, and amplitude formalism may not always be identical for different 
experiments, fit fractions are reported in addition to amplitudes to allow for more mean- 
ingful comparisons between results. The fit fraction is defined as the integral of a single 
component divided by the coherent sum of all components. The sum of the fit fractions for 
all components will in general not be unity because of the effect of interference. 

Backgrounds, combinatorics and mistags, are considered as described above. They do 
not interfere with the signal, but the mistag background shape depends on the signal shape 
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as noted earlier. 

One must also consider the statistical errors on the fit fractions. We have chosen to 
use the full covariance matrix from the fits to determine the errors on fit fractions so that 
the assigned errors will properly include the correlated components of the errors on the 
amplitudes and phases. After each fit, the covariance matrix and final parameter values are 
used to generate 500 sample parameter sets. For each set, the fit fractions are calculated and 
recorded in histograms. Each histogram is fit with a single Gaussian to extract its width, 
which is used as a measure of the statistical error on the fit fraction. 

We perform an initial fit, using the unbinned maximum likelihood technique including 
the resonances observed by E687 [11]. We then consider each of the intermediates states 
listed above retaining those that are more than three standard deviations significant. We do 
not find the k(800)-tt + , K*(U10)-it+, K 3 *(1780)-7r+ K° s p(U50), K° s p(1700), i^/ (1500), 
^5/0(1710), and the "wrong sign" _K"q(1430) + 7t~ to be significant. The Kg<j(500) is a special 
case. It is excluded in our standard fit, and its possible contribution is discussed further 
below. The remaining ten resonances, a non-resonant contribution, and backgrounds as 
described above are included in our standard fit and give our central results. 

Table I gives the results of our standard fit. Figure 2 shows the three projections of the 
fit. We note that there is a significant "wrong sign" D° — > K* + n~ amplitude. Mistags are 
not significant, having a rate of 0.1 ± 0.4%. When we compare the likelihood of our standard 
fit to one where the K* + n~ amplitude fixed to zero we see that the statistical significance 
of the K* + n~ amplitude is 5.5 standard deviations. Also we note that the phase difference 
between the K*~ir + and K* + ir~ contributions is consistent with 180°, as expected from 
Cabibbo factors. 

We consider systematic uncertainties that arise from our model of the background, the 
efficiency, and biases due to experimental resolution. Our general procedure is to change 
some aspect of our standard fit and interpret the change in the values of the amplitudes and 
phases as an estimate of the systematic uncertainty. The background is modeled with a two 
dimensional cubic polynomial and the covariance matrix of the polynomial coefficents deter- 
mined from a sideband. Our standard fit fixes the coefficients of the background polynomial, 
and to estimate the systematic uncertainty on this background shape we perform a fit with 
the coefficients allowed to float constrained by the covariance matrix. Similarly we perform 
a fit with a uniform efficiency rather than the efficiency shape determined from the simula- 
tion as an estimate of the systematic uncertainty due to the efficiency We change selection 
criteria in the analysis to test whether our simulation properly models the efficiency These 
variations to the standard fit are the largest contribution to our experimental systematic 
errors. To study the effect of the finite resolution our experiment has on the variables in the 
Dalitz plots we vary the size of the bins used to compute the overall normalization. 

Another class of systematic uncertainties arise from our choices for the decay model for 
D° — > Kgir + iT~. We consider the Zemach formalism [21], rather than the standard helicity 
model, which enforces the transversality of intermediate resonances and we vary the radius 
parameter [22] for the intermediate resonances and for the D° between zero and twice their 
standard value of 0.3 fm and 1 fm, respectively. These variations to the standard fit are the 
largest contribution to our modeling systematic errors. Additionally, we allow the masses 
and widths for the intermediate resonances to vary within their known errors [18-20]. 

We also consider uncertainty arising from which resonances we choose to include in our 
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fit to the Dalitz plot. We compared the result of our standard fit to a series of fits where each 
of the possible resonances were included one at a time. We also considered a fit including all 
possible resonances. We take the maximum variation of the amplitudes and phases from the 
standard result compared to the results in this series of fits as a measure of the uncertainty 
due to our choice of included resonances. 

The er(500) has been reported by E791 in D + — > 7r~7r + 7r + decays [23]. The parameters 
of the <t(500) are sensitive to the choice of decay model, discussed previously. Replacing 
the non-resonant contribution in our standard fit with a K° s a(500) component yields an 
amplitude of 0.57 ± 0.13 and phase of 214 ±11 with a mass of m a = 513 ± 32 and width 
T a = 335 ± 67 MeV for the a(500) consistent with E791 results [23] (m a = 478 ± 29 and 
T a = 324 ±46 MeV). While we find this suggestive that there is a K^a (500) contribution, we 
are unable to definitively confirm this because of the known shortcomings of our description 
of the scalar resonances. The systematic uncertainty does include the difference between the 
standard fit which does not include the <r(500), and the fits allowing it. On the other hand, 
we find no evidence for a scalar k~ — > Kgir~ as also suggested by E791 [24] in charm decays. 

We also do separate standard fits for D° and D° tags to search for CP violating effects. 
We see no statistically significant difference between these two fits. A more general study 
would consider a CP violating amplitude for each component observed in our standard fit. 

In conclusion, we have analyzed the resonant substructure of the decay D° — > K^it + it~ 
using the Dalitz technique. We observe ten contributions including a "wrong sign" D° — > 
K* + ir~ amplitude with a significance of 5.5 standard deviations. This decay arises from a 
double Cabibbo suppressed decay or D° — D° mixing. We measure fmoz^^+] = (0-5 ± 
0-2tai -o.i)%, and the relative phase between the two decays to be (189 ± 10 ± 3tf)°, 
consistent with 180°. We consider D° and D° tags separately, and see no CP violating 
effects. 

The time dependence of the Dalitz distribution of this decay mode in our data is able to 
discern the source of the "wrong sign" signal, doubly Cabibbo suppressed decays or D° — D° 
mixing, and have sensitivity to the mixing parameters x, explicitly to its sign, and y at the 
few percent level [25]. 

ACKNOWLEDGMENT 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 



8 



REFERENCES 



[1] M. Bauer et al, Z. Phys. C 34, 103 (1987). 

[2] P. Bedaque et al, Phys. Rev. D 49, 269 (1994). 

[3] L.-L. Chau and H.-Y. Cheng, Phys. Rev. D 36, 137 (1987). 

[4] K. Terasaki, Int. J. Mod. Phys. A 10, 3207 (1995). 

[5] F. Buccella et al, Phys. Lett. B 379, 249 (1996). 

[6] E. Golowich and A. Petrov, Phys.Lett. B 427, 172-178 (1998). 

[7] J. Adler et al, Phys. Lett. B 196, 107 (1987). 

[8] PL. Frabetti et al, Phys. Lett. B 286, 195 (1992). 

[9] J.C. Anjos et al, Phys. Rev. D 48 56 (1993). 

[10] H. Albecht et al, Phys. Lett. B 308, 435 (1993). 

[11] PL. Frabetti et al, Phys. Lett. B 331, 217 (1994). 

[12] Y. Kubota et al, Nucl. Instrum. Methods Phys. Res. A 320, 66 (1992). 

[13] T.S. Hill, Nucl. Instrum. Methods Phys. Res. A 418, 32 (1998). 

[14] R. Godang et al, Phys. Rev. Lett. 84, 5038 (2000). 

[15] G. Bonvicini et al, Phys. Rev. Lett. 82, 4586 (1999). 

[16] GEANT manual, CERN Program Library Long Writeup W5013, Copyright CERN, 

Geneva, 1993. 

[17] S. Kopp et al, Phys. Rev. D 63, 092001 (2001). 

[18] D.E. Groom et al, Eur. Phys. J. C 15, 1 (2000). 

[19] E. M. Aitala et al, Phys. Rev. Lett. 86, 765 (2001). 

[20] A. Kirk, hep-ph/0009168. 

[21] C. Zemach, Phys. Rev. B 133, 1201 (1964). 

[22] J. Blatt and v. Weisskopf, Theoretical Nuclear Physics, New York: John Wiley & Sons 
(1952). 

[23] E. M. Aitala et al, Phys. Rev. Lett. 86, 770 (2001). 

[24] C. Gobel, hep-ex/00 12009. 

[25] D. Asner, UCRL-JC-147296-DR. 



9 



600 



.400 - 



O 

> 

CD 

o 

o 
d 



1 200 

> 
LU 



t — i — i — r 



t — i — i — r 



t — i — i — r 



(a) 




1 2 

M(KgV) 2 (GeV 2 /c 4 ) 




*3280602-002 

i — i — i — i — i — i — i — i — i — 



(b) - 



J I I I I I I I L 



J I I L 



100 - 



o 

CM 
> 
CD 
O 

o 



CD 
> 
LU 



50 




1 2 ■ 

M{n + jt~f (GeV 2 /c 4 ) 

D — Jt , D — Jt Jt , Kg— +~Jt Jt 



T 1 1 1 1 1 1 1 1 1 1 1 1 T 

(d) 



6«-» i? 



J I I L 



J I I L 



J I I L 



M(KgV) 2 (GeV 2 /c 4 ) 



1 2 

(M RS ) 2 (GeV 2 /c 4 ) 



FIG. 2. Projections of the results of the fit described in the text to the K s tt + it~ Dalitz distri- 
bution showing both the fit (histogram) and the data (points). In c) the result of a fit where the 
"wrong sign" D° — ► K* + ir~ amplitude is fixed to zero is also shown, d) The Dalitz distribution 
for D° — > Kgir + ir~ candidates. The horizontal axis (Mrs) 2 corresponds to (M K o n -) for D° and 
(M K>+ )iorW. 
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TABLE I. Standard fit results. The errors shown are statistical, experimental systematic, and modeling systematic respectively. See 
the text for further discussion. 



Component 


Amplitude 


Phase 


Fit Fraction (%) 


K*(892)+vr- xB(K*(892) + -► 






(11 ±2 It -t) x W- 2 


321 ± 10 ± 3 If 


n 34 ±013 +a3i +a2y 


K°p° 






1.0 (fixed) 


(fixed) 


26.4 ±0.9 1°;? 1°; 4 


K°lo xB(lo — > 7T + 7T~) 


K°TT-) 




(37±5±1 If) x 10~ 3 


114 ± 7 ±t 11 


U.U± U.18 _ 06 _ 07 


K*(892)^vr+ xS(i^*(892)" -► 




1.56 ± 0.03 ± 0.02 tom 


150 ± 2 ± 2 11 


ar 7 4. 1 +1.1 +1.4 
DO.C zt 1.3 _ 2 6 _ 3 


K°/o(980) xS(/ (980) -► vr+vr 






34 + n 02 +0 - 04 +0 - 04 

U.34 it u.uz _ 03 _ 02 


188 ± 4 ±| ±| 


4.3 ±0.5 loi±0.5 


K°/ 2 (1270) xB(/ 2 (1270) - vr- 


f 7T-) 




0.7 ±0.2 +^±0.4 


308 ± 12 +^ +f 


27 + 1 5 + - 24 +0 - 28 

U.Zi zt U.IO _ og _ 014 


K°/o(1370) xB(/ (1370) - tt" 


f 7T-) 




1.8 ±0.118;? 


85 ±4 11 I? 4 


4- 1 1 +2-4 +1.4 
y.y zt 1.1 _ 4 3 


K,J(1430)-7r+ x 5 (if* (1430) - 


— ► A 7T 


") 


2.o±o.ii°; 2 ±8:f 


3 ± 4 ± 4l[ 5 


( .0 zt u. / _ g _ 7 


K 2 *(1430)-vr+ x5(iq(1430)- 


T7° - 
-»ft f 




1.0 ±0.1 ±0.1 ±g;f 


335 ± 7 I4 1 2 4 


11+02 +0 - 3 +a6 
1.1 zt U.Z _ ! _ .3 


K*(1680)-vr+ xB(K*(1680)- 


- 

-tft f 


") 


5.6 ±0.6 lal±4.0 


174 ± 6 ±|° +11 


2 2 + 4 + - 5 + L7 

Z.Z zt U.4 _ 3 _! 5 


-77O + - 

K 7r^7r non-resonant 






1 1 + n q +0.5 +0.9 
1.1 zt U.O _ 2 _o.7 


340 ± 11 IfJ -i 


0.9 ±0.4 1 ;° l i 



